Fast Ground-State Reordering from Off-Diagonal Fluctuations of Interaction 
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It is shown that off-diagonal fluctuations of interaction lead to the breakdown of perturbation theory 
in the immediate vicinity of the ground-state of a fermionic system at a rate inversely proportional 
to the number of considered particles. This rate is parametrically much faster than previously 
expected. It is shown that this follows from the coherent addition of many small second order 
contributions to the energy of low-lying levels. 
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Random interactions naturally occur in various fields 
of physics. Assuming the symmetry of the probability 
distribution of the one-body Hamiltonian under unitary 
basis transformation in Hilbert space directly results in 
the randomness and a similar symmetry for the distri- 
bution of the Interaction Matrix Elements (IME). Typi- 
cally, this symmetry is satisfied up to an energy scale E c 
(the Thouless energy) so that models with random inter- 
actions may be expected to describe interacting systems 
provided one restricts oneself to energy scales of an or- 
der of magnitude given by E c |l],^] . In condensed-matter 
systems, this corresponds to one-body excitation ener- 
gies inside a window of width E c = gA around the Fermi 
energy, hence containing a finite number of levels given 
by the conductance g (A is the mean level spacing). 

Random Interaction Models (RIM) first appeared in 
nuclear physics and early investigations focused on the 
deviations from Random Matrix Theory introduced by 
considering interactions of finite order k <C n (n is the 
number of particles, for k = n one recovers the Dyson en- 
sembles of random matrices) Q. The complexity of the 
excited spectrum of rare-earth atoms motivated the in- 
troduction of RIM in atomic physics Q] , and RIM repro- 
duced statistical features obtained in nuclear shell model 
calculations fl. Last but not least, RIM were used in 
condensed matter physics to study the quasiparticle life- 
time [^,0 and the conductance peak spacings and heights 
distributions H in quantum dots in the Coulomb block- 
ade regime. More recently a lot of attention has been de- 
voted to the properties of the ground-state of RIM |] Q 
and several spectral features presented in these latter in- 
vestigations indicate that RIM are realistic down to very 
low excitation energies. It is the purpose of the present 
letter to investigate in detail the ground-state structure 
of systems of few (n) fermions as the relative strength 
U/A of fluctuations of interaction increases. The main 
result of these investigations is a breakdown of pertur- 
bation theory for the ground-state at a critical value 
U c (n) ~ A/n. In particular, U c is parametrically smaller 
by a factor 1/n than previously suggested . Once this 
threshold is reached, the ground-state starts to be com- 
posed of an exponentially large number of Slater deter- 



minants, an occurrence which goes far beyond a simple 
renumbering of one-particle orbitals. We will see that 
this is due to the coherent addition of large numbers of 
small second order contributions for perturbation theory 
in the immediate vicinity of the ground-state. 
We investigate a model for spin- 1/2 fermions 
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The greek indices label m/2 different one-body ener- 
gies which are distributed as e Q S [0;m/2] so that the 
mean spacing between spin-degenerate levels is A = 1 
and ) =t, J, are spin indices. The Hamiltonian is Spin 
Rotational Symmetric (SRS): the interaction commutes 
with both the total spin \S\ and its projection S z . The 
corresponding eigenvalues a and <j z are good quantum 
numbers and the Hamiltonian matrix acquires a block 
structure where blocks are labelled by a z and a. From 
now on we consider only the a = a z = (1/2) block 
for n even (odd). The results and arguments presented 
below can however be generalized to higher spin blocks 
or models of spinless fermions. 

The Hamiltonian (1) can be viewed as a model of inter- 
acting fermions expressed in the basis of Slater determi- 
nants constructed from the set of one-particle eigenfunc- 



tions 1p a of H = e aCa, 

are given by U^'g 



,. In this basis, the IME 

where U(i — j) is the interaction potential in real space; 
i and j label sites on a lattice. Assuming that the pres- 
ence of disorder or chaotic boundary scattering confers 
a random character to the ^> Q 's directly leads to fluc- 
tuations of the U^'p around their average value. Only 
diagonal matrix elements have a nonzero average leading 
to mean-field charge-charge, spin-spin and BCS interac- 
tions jlo| which we will neglect as they have at most a 
marginal influence on the object of interest here. We 
therefore take the IME in (1) to be randomly distributed 
inside a zero-centred Gaussian distribution of width [/, 
P(l^)oce^)V^\ 

We start with a perturbative calculation of the ground- 
state energy up to the second order in the small param- 
eter U/A 1. The first order correction vanishes on 
average, and the second order correction can be written 
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e<5 



e/3 



(2) 



indicates that the sum is taken over occupied lev- 
els 7 and S and unoccupied ones a and (3 (taking into 
account spin degeneracy), giving a total number of con- 
tributions K ~ n 2 (m, — n) 2 . For the ground-state, one 
has e 7 + e<5 < e Q + cp, and each term in the sum (||) is 
negative. The contributions from each scattering process 
therefore add coherently resulting in a strong reduction 
of the ground-state energy. The latter can be estimated 
for an equidistant spectrum by replacing sums by inte- 
grals in (Q). This gives an homogeneous polynomial of 
order three in n and m, each term being multiplied by a 
logarithmic correction. In the dilute limit 1 <C n <C m 
the m 3 and m?n terms drops out exactly so that the 
dominant contribution to the second order correction in 
the ground-state energy reads 



AS, 



(2) 



-A{U 2 /A)n 2 mlog(m) 



(3) 



The prefactor A can be estimated A r* 7; details of the 
perturbative treatment are presented elsewhere p0[. 
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FIG. 1. Average ground-state energies for 500 realizations 
of Hamiltonian (1) with m = 16 orbitals, n = 2 (circles), 4 
(squares), 6 (diamonds) and 8 (triangles) compared to the sec- 
ond order perturbative result (3) (dashed lines) and to linear 
fits valid at large U (solid lines - see text). 

In Fig.l we show the behavior of the ground-state en- 
ergy for n — 2, 4, 6 and 8. It is clearly seen that the av- 
erage ground-state energy exhibits the quadratic depen- 
dence in U predicted by (2-3) at weak U, while at larger 
U, this dependence becomes linear. Most importantly, 
this crossover occurs at a weaker U for larger n. Both 
these features can be understood as one realizes that for 
dominant fluctuations of interaction U/A ^> 1 the Hamil- 
tonian depends linearly in U. In this regime, and for 
n ^> 1, the average many-body density of states is Gaus- 
sian with a width B well approximated by B ~ nraU 



0,0. This width thus gives a good parametric esti- 
mate for the average ground-state energy in the asymp- 
totic regime ^q 2 "* ~ B ~ nmU || so that equating the 
asymptotic and perturbative expressions for the ground- 



state energy (A£q 



(2) 



«j nmU) gives a critical value 
U c ~ A/n (4) 



The crossover between the two regimes occurs at a rate 
given by U c and is inversely proportional to the number 
of particles. As can be seen in Fig.l, below this thresh- 
old, the ground-state energy is well approximated by the 
second order perturbation expression (3), while above it, 
its [/-dependence becomes linear which cannot be cap- 
tured by any finite order of perturbation theory, and we 
are forced to conclude that U c gives the radius of conver- 
gence for perturbation theory in U. 

For U > U c , the ground-state starts to be mixed 
with higher excited states and the number of Slater 
determinants J) necessary to construct it becomes 
much larger than one. A measure of this number 
is given by the Participation Ratio (PR), which, for 
a given many-body eigenstate \^\(U)), is defined as 

& = (E,7I^I*a(c/)}| 2 ) 2 /EjI<^I*a(^)}| 4 . if \9x(U)) 

is made up of one single Slater determinant (as it is for 
U = 0), one has £\ = 1 whereas when \ty\(U)) is homo- 
geneously spread over the N Hilbert space basis states, 
£a ~ N. To check the above threshold (0), we numer- 
ically calculated the PR for the ground-state wavefunc- 
tion \^>q(U)). Results obtained from the exact diagonal- 
ization of small systems of n < 8 particles are shown in 
Fig. 2. We first note a faster increase of £o with U at 
larger n which is not due to an increase of the Hilbert 
space size as it occurs before a significant portion of the 
Slater determinant basis starts to be occupied. Most im- 
portantly, this increase is exponential (note the vertical 
log-scale on Fig. 2) : fixing e.g. U/A = 0.2, one gets 
£(n) w exp(n/1.6). Since each order in perturbation the- 
ory gives only an algebraic number of contributions, this 
exponential increase of the PR cannot be captured by 
a finite, n-independent number of orders, and therefore 
indicates the divergence of the perturbation series as n 
increases. Another important feature of Fig. 2 is the quite 
small value U <Si A at which £o starts to increase. This 
is quite surprising as A is the only energy scale close to 
the ground-state. When the PR increases beyond £o > 2, 
second (and higher) order terms become more and more 
important, until eventually perturbation theory breaks 
down. One can thus extract parametrically the radius 
of convergence of perturbation theory from the PR by 
defining a critical interaction strength U c from the con- 
dition £o(U c ) = 2. (This choice is arbitrary as long as 
one chooses a threshold value £o(U c ) much smaller than 
the total size of the Hilbert space as this would intro- 
duce an artificial n and m dependence in U c .) The inset 
to Fig. 2 shows a n— (and absence of m— ) dependence of 
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U c clearly confirming 
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FIG. 2. Ground-state PR £ vs. U/A averaged over 100 
(for n = 8) to 2000 (for n = 2) realizations of Hamiltonian 
(1) for m = 16 and n=2, 3, ... 8 (solid lines from bottom to 
top). The dotted dashed-line indicates the threshold £o = 2. 
Inset : Critical value of interaction U c /A necessary to reach 
£o = 2 vs. n, for m=10 (triangles), 12 (+), 14 (x) and 16 
(circles). The dashed line gives a decay U c /A = 0.35/n. 

As yet another check of the threshold (f|), we consider 
the onset of level repulsion for the spacing si = £\ — £q 
between the ground-state energy £q and the first ex- 
cited level £\. At U = 0, this spacing corresponds to 
a one-body excitation, and considering a randomly dis- 
tributed one-body spectrum, it is Poisson distributed 
Pp{ s i) = ex p(— Switching on U induces off-diagonal 
matrix elements connecting the two lowest levels and 
leads to repulsion between them. Making a two-level ap- 
proximation incorrectly gives the onset of this repulsion 
at U = A. A correct treatment goes as follows. Most 
terms in the second order contribution (|^) are common to 
both the energies of the ground- and first excited states. 
The difference corresponds to scattering from or onto the 
(n/2) th and (n/2 + l) th orbitals, giving a contribution to 
the energy difference between the lowest two levels 
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/2+l> 



£6 
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(There is an additional similar contribution arising from 
the scattering from the n/2 and n/2+1 orbitals.) We first 
make the substitution a((U^ n/2 ) 2 - (U^ +1 ) 2 ) = U 2 
for IME fluctuations. Then, and considering an equidis- 
tant spacing, a simple counting shows that the number of 
contributions in the above expression corresponding to a 
given energy denominator pA grows as p 2 up to p = n/2 
and decreases afterward. Neglecting logarithmic correc- 
tions, one gets an estimate for the contributions to the 
relative fluctuations given by - (U 2 /A)n 2 . This estimate 
can also be obtained by replacing the sums in (§) with 



integrals. Once these fluctuations become of the order 
of the spacing between £ and £%, — (U 2 /A)n 2 — A, an 
avoided crossing occurs. This signals the onset of level 
repulsion which therefore occurs at the rate given by (H). 
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FIG. 3. (a,b,c) Distribution -P(si) of the first level spacing 
for 50000 realizations of Hamiltonian (1) with a random dis- 
tribution of orbitals e a , m = 14 (symbols) n = 3, 4, ... 7, at 
an interaction strength U/A = 1/n. (a) Shows P(si) on a 
normal scale, (b) its large si tail on a semilogarithmic scale 
and (c) its small si behavior on a log-log scale. Empty sym- 
bols correspond to odd n and full symbols to even n. Clearly, 
there is an even-odd dependence which is m-independent as 
shown by the superimposed values for m = 18 and n = 3 
(dashed line) and n = 4 (dotted-dashed line) on (b) and (c). 
On (a,b,c) the solid lines show the Poisson and Wigner-Dyson 
distributions, (d) Shows one realization of the spectrum for 
Hamiltonian (1) with equidistant one-body spectrum. The 
spectrum correspond to the full sector a z = so that level 
crossings occur between levels with different a. 

This reasoning is confirmed on Fig. 3 (a,b,c) where we 
present the distribution P(si) of the first level spacing 
for Hamiltonian (1) and different number of particles n 
at fixed interaction strength U = A/n. The two main 
features that can be seen are : 1) P{s±) at fixed U = U c 
is n-independent, except that 2) P(si) shows a clear 
even-odd dependence, in particular there is more level 
repulsion for odd n. This second fact is understood once 
one realizes that due to SRS, and within a given (er 2 , a)- 
block, the first U = excitation is doubly degenerate for 
n odd, while it is not degenerate for n even. This results 
in a stronger effect of the interaction in the former case 
and the even-odd effect seen on Fig. 3 (a,b,c). It is to 
be stressed that this is purely a spin effect, this even- 
odd dependence does not appear for spinless fermions. 
On Fig. 3 (d) we also show a plot of the low-lying part 
of the spectrum which shows an avoided crossing at a 
value U/A w 0.05 much too small to be explained by 
a two-level approximation. We stress that in (d) the 
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avoided crossing is not accompanied by an increase of 
the level spacing as would be the case in presence of non- 
fluctuating IME, and which would delay the breakdown 
of perturbation theory. 

Equation (Q) predicts a parametrically small radius of 
convergence for perturbation theory beyond which the 
ground-state is composed of exponentially many Slater 
determinants. It is somehow in opposition with earlier 
works which, based on extrapolations from higher 
excitation energies, predicted a border ~ A. The mecha- 
nism leading to the breakdown of perturbation theory at 
higher excitation energies is however due to mixing be- 
tween two-body levels which are very close in energy, and 
thus give the dominant contribution in a second order 
approach. In the immediate vicinity of the ground-state 
these processes with small denominators do not exist as 
the lowest energy scale is A - the approaches used in Q 
correspond then to a two-level approximation involving 
ground- and first excited states. The presented theory 
show that the correct procedure is to retain a large num- 
ber of small second order contributions which have an 
enhanced effect close to the ground-state, as they acquire 
a well-defined sign and hence add coherently to give the 
small threshold (pj) . One has to conclude that the ap- 
proach used in |12[| , while fully justified at a sufficient 
excitation energy, must break down somewhere close to 
the ground-state. 

Random interaction models of the kind investigated 
here describe interacting systems restricted to energy 
differences of the order of the Thouless energy E c Q. 
In disordered diffusive systems one has E c = gA (g is 
the conductance) and introducing a cut-off at E c leads 
to to, n ~ g. From (Q) one would expect the break- 
down of perturbation theory to occur for U c ~ A/g, pre- 
cisely at the amplitude of fluctuations of IME in these 
systems |]J^|. The studied model neglects however 
two ingredients which must be considered when deal- 
ing with more realistic systems. First, correlations be- 
tween the two terms in (^J) significantly reduce the value 
of this latter expression below the uncorrelated estimate 
U 2 ~ (A/g) 2 . Second, the mean-field charge-charge in- 
teraction (that we voluntarily neglected), when incorpo- 
rated self-consistently into the one-body spectrum, leads 
to a strong reduction of the one-body density of states 
fL6f , hence an increase of the level spacing at the Fermi 
level against which (^) must be weighted, and an increase 
of the denominators for the infrared contributions to sec- 
ond order perturbation theory. It is expectable that these 
two facts ensure the convergence of perturbation theory 
and the stability of the ground-state against off-diagonal 
fluctuations in diffusive metallic systems. 

In summary, based on second order perturbation the- 
ory we predicted a fast reordering of the ground-state for 
systems of randomly interacting fermions which we con- 
firmed numerically. While the influence of this reorder- 
ing on the quasiparticle excitations is at present unclear, 



we found that it is accompanied by a change of statis- 
tics of the first energy excitation above the ground-state. 
Similar investigations in lattice models where this onset 
of level repulsion may influence the system's transport 
properties pTj , and for which correlations between IME 
should play an important role (neglected in the presently 
studied model) are highly desirable. 

Numerical computations were performed at the Swiss 
Center for Scientific Computing. Work supported by the 
Swiss National Science Foundation. 
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